The experimental observation of Beliaev damping in a Bose condensed gas 
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We report the first experimental observation of Beliaev 
damping of a collective excitation in a Bose-condensed gas. 
Beliaev damping is not predicted by the Gross-Pitaevskii 
equation and so this is one of the few experiments that tests 
BEC theory beyond the mean field approximation. Measure- 
ments of the amplitude of a high frequency scissors mode, 
show that the Beliaev process transfers energy to a lower lying 
mode and then back and forth between these modes. These 
characteristics are quite distinct from those of Landau damp- 
ing, which leads to a monotonic decrease in amplitude. To 
enhance the Beliaev process we adjusted the geometry of the 
magnetic trapping potential to give a frequency ratio of 2 to 
1 between two of the scissors modes of the condensate. The 
ratios of the trap oscillation frequencies uj v /uj x and lj z /uj x 
were changed independently, so that we could investigate the 
resonant coupling over a range of conditions. 

PACS numbers: 03.75.Fi, 05.30.Jp, 05.45.-a, 32.80.Pj 



In the Beliaev damping process [Q, one quantum 
of excitation converts into two quanta of a lower fre- 
quency mode, in a manner analogous to parametric 
down-conversion of light in nonlinear media. This process 
occurs readily in a homogeneous system such as super- 
fluid helium where there are many very closely spaced 
low- lying energy levels. However, it has not previously 
been observed for a trapped Bose-condensed gas, which 
has discrete modes with well-resolved excitation energies. 
We have observed the Beliaev process for a scissors mode 
H of a BEC of rubidium atoms, when the mode res- 
onantly couples to another mode at half its frequency. 
Beliaev damping can be regarded as a four wave mix- 
ing process, in which a Bogoliubov quasiparticle inter- 
acts with the ground state component of the condensate 
to produce two new quasiparticles that share the initial 
energy equally. 

At zero temperature, only Beliaev damping may occur, 
whilst at finite temperature Landau damping is the pri- 
mary mechanism which dissipates the energy of collective 
excitation into the thermal cloud [Q. Landau damping 
corresponds to a scattering process between two quasi- 
particles in the initial mode, which results in one par- 
ticle gaining all of the energy and the other going into 
the condensate ground state. The Landau damping rate 
and its strong dependence on temperature has been mea- 
sured for two of the quadrupole excitations of the con- 
densate We observe a comparable Landau damping 



rate for the scissors mode in most trap geometries, ex- 
cept those where the frequencies of the two scissors modes 
were matched and the Beliaev process was dominant. 

An atom in a harmonic trap has three dipole modes of 
oscillation along each of the principle axes, corresponding 
to the three trap frequencies u) x , uj y and u> z . The centre of 
mass motion of an interacting cloud of atoms, such as the 
Bose condensate and any surrounding thermal cloud, has 
the same eigenfrequencies by Kohn's theorem ||. How- 
ever the frequencies of other collective excitations of the 
condensate, including the quadrupole modes, do not cor- 
respond to those of a harmonic oscillator because of the 
strong interatomic interactions. Six of the low lying col- 
lective excitations are described by the linearised set of 
equations: 



q = Sq 

where the six components of the vector qf, 

Q {Qiii Q jj ? Qkki Qij i Q jki Qik) 

are elements of the quadrupole tensor, Q. 
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The six eigenvalues of S give the frequencies of the 
six normal modes. Linear combinations of the diagonal 
elements of Q ( e.g. (x 2 )) describe the three normal 
modes with principal axes of fixed orientation. We will 
refer to these modes as Mh,,Mi,Mz. In the case of the 
axially symmetric trap, where the angular momentum 
about the axis (m) is a good quantum number, these 
become the high and low lying m = modes and the 
m = 2 mode respectively The three independent 

off-diagonal elements of Q (e.g. (xy)) relate to the three 
scissors modes, which we label M xy , M yz , M xz , with fre- 
quencies uj xy , ujy Z , uj xz . The frequency spectrum has been 
calculated in and is shown, in the hydrodynamic 

limit, in Fig. pTas a function of trap geometry. 

In general, the first and second harmonic oscillations 
involved in the down-conversion process are the scissors 
modes in the xy and xz planes respectively (M xy and 
M xz ). However, for the special case of an axially sym- 
metric trap, M xy is degenerate with Mi and the two 
cannot be distinguished. The frequencies of these two 
scissors modes in the hydrodynamic limit are given by 
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ui xy = y^x + and uj xz — \/w x + w z . Thus we ob- 
serve a resonant coupling when the following condition is 
satisfied: 

The scissors mode is a small angle, irrotational oscil- 
lation of the condensate, about a given axis, that occurs 
without a change of the cloud shape p|JTl|] . Mi is an out 
of phase oscillation in the x and y directions, with a very 
small motion in the z direction. This axial motion falls to 
zero amplitude in the cylindrically symmetric trap. Ref- 
erences show the geometries of the two oscillations 
& 

We initially create Bose-Einstein condensates in an 
axially symmetric time-orbiting potential (TOP) trap 
(lu x = ujy = 127.6 ± 1.0 Hz). The trap frequency ratio, 
uj z /uj x , is 2.83 and hence the condensate has an oblate 
shape. After RF evaporative cooling, the condensate 
temperature is significantly below 0.5 T c and the ther- 
mal cloud is no longer visible. A detailed description of 
our method for selectively exciting the xz scissors mode 
is given in jll[] . In summary, we apply an additional bias 
field of amplitude B z in the axial direction, oscillating 
in phase with the x component of the TOP bias field, 
B x . This enables the symmetry axis of the trapping po- 
tential to be tilted by an angle that depends on B z /B x . 
We first form the condensate in a standard TOP trap, 
then adiabatically tilt the trap to an angle <j> and then 
suddenly flip it to — (j). This excites a scissors mode os- 
cillation of amplitude 9 — 2(f> in the xz plane, about the 
new equilibrium position. 

However the application of the z bias field has a second, 
more subtle effect The axial trapping frequency, 

uj z , and hence the xz scissors mode frequency is reduced, 
as B z is increased. The radial frequencies u) x and u> v are 
also reduced, but by a negligable amount. By changing 
B z j B x we are able to tune the xz scissors frequency into 
resonance with twice the xy scissors frequency. Since the 
excitation angle is coupled to B z /B x , this is also changed, 
but always remains well within the small angle scissors 
mode limit defined in j|. Figure || shows theoretically 
and experimentally how the excitation amplitude, 9, and 
the frequency of the xz scissors mode, lo xz , change as a 
function of B z /B x and hence uj z /uj x . 

Our first observation of Beliaev damping occurred in 
an axially symmetric trap, tu x = u) y . In this case the 
xy scissors mode becomes degenerate with, and indistin- 
guishable from the \m = 2| quadrupole mode, with fre- 
quency y/2 u> x (Fig. [I]) . The resonance condition of Eq. 
([|) is then fulfilled when lo z /uj x — \fl. 

We recorded the xz scissors mode oscillation using de- 
structive absorption imaging. The angle of the cloud as a 
function of oscillation evolution time, was extracted from 
2D gaussian fits to absorption images of the expanded 



condensate. We fit an exponentially decaying sine wave 
to the angle versus time data and use the damping rate 
as a measure of the rate at which energy is flowing out 
of the mode. For values of lo z /uj x far from the resonance 
condition, we observe only a slow Landau damping rate 
of ~18 Hz (Fig. ||a). This background damping rate 
increases with the mode frequency, w I2 , as predicted in 

Close to resonance we observe two new features in the 
xz scissors mode. Firstly there is a sharp increase in the 
initial damping rate of the oscillation, as the Beliaev pro- 
cess becomes significant (Fig. ||b). Secondly, over longer 
observation periods ( t > 30ms ), the amplitude enve- 
lope grows and then falls again, as energy continues to 
be transferred back and forward between the two modes, 
in a manner characteristic of non-linear coupling (Fig. 
|c) To obtain a measure of the total damping rate, 
we fit the decaying sine wave up to, but not beyond, the 
first amplitude minimum. The maximum coupling rate 
occurs when uj z /uj x = 2.65±0.04, in very good agreement 
with the predicted value of \fl (Fig. ||). 

In the axially symmetric trap described so far, Mi and 
M xy are degenerate, but in a totally anisotropic trapping 
potential, lu x ^ uj y , this degeneracy is broken (Fig. [l]). 
To investigate coupling to the pure xy scissors mode, we 
repeated the experiment, looking for resonant behaviour 
in a range of triaxial traps. The trap geometry was cho- 
sen to be close to the resonance condition between M xy 
and M xz given in Eq. (j|). 

The condensate was formed in a symmetric TOP trap 
as before. However, whilst it was slowly tilted in the 
xz plane, it was also adiabatically deformed by ramp- 
ing the amplitude of the y component of the TOP bias 
field. This made the trapping potential elliptical in the 
xy plane, with uj y /uj x between 0.6 and 1.05 fl6T| . For 
each value of uj y /oj x , we determined the value of lo z /lo x 
for peak coupling. Figure [| shows this data, with the res- 
onant coupling condition superimposed. The agreement 
between experiment and theory over a range of trap ge- 
ometries confirms that the correct damping process has 
been identified. 

These observations do not rule out the possibility of 
down-conversion into the Mi mode in an anisotropic po- 
tential, since our traps were specifically chosen for M xz 
to be resonant with M xy and not Mi. In fact, for traps 
with a very small anisotropy in the xy plane, so that 
both modes are close to resonance with the xz scissors 
mode, we notice a broadening of the peak in the initial 
damping rate. This is expected when a second down- 
conversion process, into the Mi mode, is resonant in a 
slightly different trap geometry. 

In non- linear optics, it is well-known that down- 
conversion produces output beams of squeezed light and 
analogous effects with matter waves may be explored us- 
ing BEC. In previous work we observed up-conversion 
between the two m — modes of a BEC in a cylindri- 
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cally symmetric potential [ |l7fl and interesting compar- 
isons may be made between the two processes. Whilst 
the up-conversion process showed strong dispersion of the 
driving oscillation close to resonance, this was not ob- 
served in the measured values of uj xz in Fig. ^. Secondly 
the return of energy to the original mode, shown strik- 
ingly in Fig. ^Jc, was not observed for the up-conversion 
process. Finally, whilst spontaneous up-conversion is 
predicted by the Gross-Pitaevskii equation, the reverse 
process of spontaneous down-conversion may not be de- 
scribed within this framework. This may be explained by 
the following simple argument |L5]]. The GPE contains 
a nonlinear term proportional to |i/ , | 2 J so excitation of a 
pure mode with time dependence e 4 "', leads to terms of 
angular frequency 2u>, but not to terms in to/ 2 and hence 
no down-conversion. However, if there is initially some 
amplitude at frequency w/2, then nonlinear mixing with 
the oscillation at frequency u> leads to a transfer of energy 
between the modes at uj and uj/2. 

In conclusion, we have observed Beliaev damping of the 
xz scissors mode. Measurements over a range of trap ge- 
ometries confirm that this corresponds to resonant down- 
conversion into the xy scissors mode. In future experi- 
ments, we intend to probe the nature of the final squeezed 
state and investigate the role of the thermal cloud in the 
down-conversion process at finite temperature. 

We would like to thank all the members of the Ox- 
ford theoretical BEC group, in particular K. Burnett, M. 
Rusch and S. Morgan for their help and advice. 

This work was supported by the EPSRC and the TMR 
program (No. ERB FMRX-CT96-0002). O.M. Marago 
acknowledges the support of a Marie Curie Fellowship, 
TMR program (No. ERB FMBI-CT98-3077). 



[1] S.T. Beliaev, Sov. Phys. JETP. 34, 323 (1958) 

[2] M.H. Anderson et a/., Science 269, 198 (1995); K.B. 

Davis et ai, Phys. Rev. Lett. 75, 3969 (1995). 
[3] D. Guery-Odelin and S. Stringari, Phys. Rev. Lett. 83, 

4452 (1999) 
[4] S.A. Morgan J.Phys.B, 33 3847 (2000) 
[5] D.S. Jin, M.R. Matthews, J.R. Ensher, C.E. Wiemann, 

E.A. Cornell, Phys. Rev. Lett. 78, 764 (1997); D.M. 

Stamper-Kurn, H-J. Meisner, S. Inouye, M.R. Andrews, 

W. Ketterle, Phys. Rev. Lett. 81, 500 (1998) 
[6] W. Kohn, Phys. Rev. 123 1242 (1961) 
[7] D.S. Jin, J.R. Ensher, M.R. Matthews, C.E. Wieman and 

E. A. Cornell, Phys. Rev. Lett. 77, 420 (1996); 
[8] M.O. Mewes, M.R. Andrews, N.J. van Druten, D.M. 

Kurn, C.G. Townsend and W. Ketterle, Phys. Rev. Lett. 

77, 988 (1996) 

[9] F. Dalfovo, S. Giorgini, L.P. Pitaevskii and S. Stringari, 
Rev. Mod. Phys. 71 463 (1999) 



[10] 
[11] 

[12] 



[13] 
[14] 

[15] 
[16] 

[17] 



M.O. da C. Pires and E.J.V. de Passos, J.Phys.B. 33 
3929 (2000) 

O.M. Marago, S.A. Hopkins, J. Arlt, E. Hodby, G. 
Hechenblaikner and C.J. Foot, Phys. Rev. Lett. 84, 2056 
(2000) 

Note that in an axially symmetric trap down-conversion 
occurs into an equal superposition of the m = +2 and 
m — — 2 modes, ensuring that the net angular momentum 
about the z axis is zero at all times. For a triaxial trap, 
angular momentum is not conserved and need not be 
considered in the coupling process. 

E. Hodby, G. Hechenblaikner, O.M. Marago, J. Arlt, S. 

Hopkins and C.J. Foot, J.Phys.B. 33, 4087 (2000) 

L.P. Pitaevskii and S. Stringari, Phys. Lett. A 235, 398 

(1997); P.O. Fedichev, G.V. Shlapnykov, and J.T.M. 

Walraven, Phys. Rev. Lett. 80, 2269 (1998) 

S.A. Morgan, S. Choi, K. Burnett and M. Edwards, Phys. 

Rev. A, 57, 3818 (1998) 

J.R. Ensher, PhD thesis, University of Colorado (1998); 
J. Arlt, O. Marago , E. Hodby, S.A. Hopkins, G. Hechen- 
blaikner, S. Webster and C.J. Foot, J. Phys. B 32, 5861 
(1999) 

G. Hechenblaikner, O.M. Marago, E. Hodby, J. Arlt, S. 
Hopkins and C.J. Foot, Phys. Rev. Lett. 85, 692 (2000) 



I 

I 

I 





_ M 

xz ... — ----- 


M 

> 


I 


< 

A 

1 1 1 


r 



1.2 1.4 



CO / CO 

y x 



FIG. 1. The frequency spectrum of the six quadrupole 
modes of the condensate as a function of the ellipticity of 
the trapping potential in the xy plane. lo z /lu x is chosen to be 
\/7, at which value M xz has exactly twice the frequency of 
M xy and M2, in an axially symmetric trap. 



3 



0.52 


B /B 

Z X 

0.79 1.08 


1.40 


- ® 








.-5- 






- 2. 






















r' 

z 






5 






- 5 






2.8 2.6 


2.4 2.2 


2.0 



2» CD 
CD 

era 

10 ^— * 



Trap frequency ratio, co z / co x 



FIG. 3. Experimental data showing the angle of a conden- 
sate in the xz scissors mode as a function of time, (a) and 
(b) are typical plots from cylindrically symmetric traps with 
different aspect ratios, fitted with exponentially decaying sine 
waves, from which the damping rates, T, were extracted. 

(a) <j z /(J x = 2.36, T = 17 ± 5 Hz. 

(b) uj z /lu x = 2.62, r = 57 ± 8 Hz. 

In (b) the trap geometry is very close to the resonance condi- 
tion of Eq. ^. and the fast decay is primarily due to Beliaev 
damping. The conditions of (a) are far from resonance and so 
only the slow underlying Landau damping rate is observed. 
Figure (c) has the same initial data as (b) but continues to 
track the oscillation over a longer period. The transfer of 
energy back and forth between the two modes is clearly ob- 
served, as the amplitude of M xz falls and rises. 



FIG. 2. The various effects of increasing the z bias field, B z . 
B z /B x is shown along the top axis and the trap frequency ra- 
tio, lu z /lj x , that results is shown along the lower x axis. The 
solid line is the theoretical value of the xz scissors frequency, 
uo xz , scaled against lu x . The solid circles show experimental 
values obtained by fitting sine waves to data such as in fig. 
[|i,b (left axis). The dotted line shows the theoretical excita- 
tion angle, 6, whilst the open circles show the values obtained 
from measuring the tilt angle of the static trap (right axis). 
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FIG. 4. Damping rate of the xz scissors mode as a function 
of the trap frequency ratio ui z /lu x , in an axially symmetric 
trap. The data is fitted with a Lorentzian peak plus a linear 
function. The Lorentzian gives the trap geometry for peak 
Beliaev damping. Far from resonance, the Landau damping 
rate is observed. The increase of the Landau damping rate 
with mode frequency is modeled by the linear part of the fit. 
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FIG. 5. Trap geometries for which M xz is resonant with 
twice the frequency of M xy (solid line) and M2 (dashed line) . 
The data points show the value of lu z /(±>x at which the peak 
damping rate was observed, for a given uj v /oj x . 



